Matter density perturbations and effective gravitational constant 
in modified gravity models of dark energy 
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We derive tlie equation of matter density perturbations on sub-liorizon scales for a general La- 
grangian density f{R, (f>, X) that is a function of a Ricci scalar R, a scalar field and a kinetic term 
X = — (V0)^/2. This is useful to constrain modified gravity dark energy models from observations 
of large-scale structure and weak lensing. We obtain the solutions for the matter perturbation 5,n 
as well as the gravitational potential $ for some analytically solvable models. In a f{R) dark energy 
model with the Lagrangian density f{R) — aR^^"^ — A, the growth rates of perturbations exhibit 
notable differences from those in the standard Einstein gravity unless m is very close to 0. In scalar- 
tensor models with the Lagrangian density / = F{(f))R-\-2p{4>,X) we relate the models with coupled 
dark energy scenarios in the Einstein frame and reproduce the equations of perturbations known 
in the current literature by making a conformal transformation. We also estimate the evolution 
of perturbations in both Jordan and Einstein frames when the energy fraction of dark energy is 
constant during the matter-dominated epoch. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 



Recent observations have determined basic cosmologi- 
cal parameters in high-precisions, but at the same time 
they posed a serious problem about the origin of dark 
energy (DE). The analysis of Super-Nova la (SNIa) [l| 
is based upon the background expansion history of the 
universe around the redshift z < 0(1). The constraint 
obtained from SNIa so far has a degeneracy in the equa- 
tion of state (EOS) of DE 0. To many people's frus- 
tration, the ACDM model with an EOS ujde = ^ 1 has 
been continuously favored from observations. This de- 
generacy has been present even adding other constraints 
coming from Cosmic Microwave Background (CMB) Q 
and Baryon Acoustic Oscillations (BAO) [J- 

The models of dark energy can be broadly classified 
into two classes [1, @ . The first corresponds to introduc- 
ing a specific matter that leads to an accelerated expan- 
sion. Most of scalar field models such as quintessence 
and k-essence [1] belong to this class. The second 
class corresponds to so-called modified gravity models 
such as /(i?) gravi ty scalar-tensor theories [l^l and 
braneworld models . In order to break the degeneracy 
of observational constraints on iude and to discriminate 
between a host of DE models, it is important to find addi- 
tional information other than the background expansion 
history of the Universe. In this paper we will show that 
modified gravity models can be distinguished from others 
by considering the evolution of matter perturbations 5m 
and gravitational potentials <f>. 

In Einstein gravity it is well known that linear matter 
perturbations on sub-horizon scales satisfy the following 



equation 



5„ 



(1) 
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where H is a. Hubble parameter, G is a Newton's grav- 
itational constant, pm is an energy density of the non- 
relativistic matter, and a dot represents a derivative with 
respect to cosmic time t. During the matter-dominated 
epoch this has a growing- mode solution oc a cx t^/^, 
which leads to the formation of large-scale structure. 
In modified gravity models the growth rates of pertur- 
bations are different because of the modification of the 
gravitational constant as well as the change of the back- 
ground evolution. In the context of f{R) gravity, in par- 
ticular, there have been a number of recent works about 
the evolution of density perturbations during the matter- 
dominated and dark energy dominated epochs [l3 |. 

We will derive the equation of matter perturbations in 
Sec. Ill for a very general Lagrangian density /(i?, (f>, X), 
where i? is a Ricci scalar and (/) is a scalar field with a ki- 
netic term X = — (V<^)^/2. Together with a sub-horizon 
approximation we assume that F = df /dR depends on 
(j) and R but not on X. In fact this Lagrangian cov- 
ers most of modified gravity DE scenarios such as f{R) 
gravity models and scalar-tensor theories. The effect of 
modified gravity appears in an effective gravitational con- 
stant Geff whose explicit form is given in Eq. (|iD|) . We 
derive a parameter rj introduced in Ref. to quantify 
the strength of an anisotropic stress and also evaluate a 
parameter E = 5(1 -I- 77/2), where g is a quantity that 
characterizes the deviation from the gravitational con- 
stant measured in solar system experiments today. 

The results in this paper can be important for future 
surveys of weak lensing as well as for the observations 
of large-scale structure (LSS) In Ref. [3 the de- 

viation from Einstein gravity was constrained from the 
galaxy clustering by taking into account an additional 
Yukawa correction to the gravitational constant. It will 
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be possible to carry out similar observational constraints 
on our f{R, (j), X) DE models from the LSS data by solv- 
ing the equation of matter perturbations. In Ref. the 
authors proposed a DE parametrization using the vari- 
ables (Sjjy) together with a linear perturbation growth 
factor 7 introduced in Refs. [13, E, If the devia- 
tion from the Einstein gravity case (S, 77) = (1,0) is de- 
tected from future survey of weak lensing, this allows us 
to distinguish modified gravity models from the models 
in Einstein gravity. 

In Sec. IV we will find solutions for 6m and $ during 
the matter-dominated epoch for some analytically solv- 
able models. In particular we show that f{R) dark energy 
models have a peculiar scale-dependence of perturbations 
unlike the case of Einstein gravity. The effect of modi- 
fied gravity on perturbations is important provided that 
a dimensionless variable m = Rf ^rrJ f .r-, which charac- 
terizes the deviation from the ACDM model, is not very 
close to 0. 

The scalar-tensor models with the Lagrangian den- 
sity / = F{(j))R + 2p{4), X) correspond to coupled dark 
energy models in the Einstein frame with a coupling 
Q{4>) = ~F_tf,/2F. In Sec. V we derive the equation of 
matter perturbations in the Einstein frame under a con- 
formal transformation and show that this in fact coin- 
cides with the equation in the models of dark energy cou- 
pled to the matter [iO]. We also derive the growth rates 
of perturbations in both Jordan and Einstein frames for 
the models in which the so-called 0-matter dominated 
epoch [m is present. 



equations 

3FH^ 



f,xX + -{FR- f)-3HF 



~2FH = f^xX + F-HF + P, 
^(a'0/,x)'-/,0 = O, 

Pm + SHpm = , 



(5) 
(6) 
(7) 
(8) 



where H = a/a, R = 6{2H^ + H), and a dot represents 
a derivative with respect to cosmic time t. 

In order to confront the DE equation of state with 
observations such as SNIa. we rewrite Eqs. ^ and ^ as 
follows: 



ZFqH — + Pm , 

-2FqH = pde +PDE 



(9) 
(10) 



where 



PDE = ^{FR - /) - iHF + f,xX + SH^iFo ~ F),(ll) 

PUE = F + 2HF -l{FR-f)- {2H + 3H^)iFo - F). 

(12) 

Here the subscript "0" represents present values. It is 
easy to show that pde and pde defined in this way satisfy 
the usual energy conservation equation 



Pde + 3i? (pde + Pde) = , 



(13) 



II. BACKGROUND EQUATIONS 



We start with the following 4-dimcnsional action 



S 



-f{RA,X)+Cr, 



(2) 



where g is a determinant of a metric 5^,^, / is a function 
in terms of a Ricci scalar i?, a scalar field and a kinetic 
term X = — ^''^i^^c/S. Cm is a Lagrangian density for a 
pressureless matter whose energy density is given by pm- 
We use the metric signature (—,+,+,+). 

The gravitational field equation and the equation of 
motion of the field (p are given by 



FG^, = \{f- RF)g^, + F,,;, - UFg^, 



+ -^.f.x(l),tj.(l),i 

(/,X0'^);c + /0=O, 



(3) 
(4) 



where we used Eq. ([7]). This was already shown to hold 
in the context of the scalar-tensor gravity [13, [1^ as well 
as the f{R) gravity [l^l- We define the DE equation of 

state as 



Pbe 
Pde 

-1 + 



2f^xX + 2F - AHF - 4H{Fo - F) 
2f^xX + FR- f- 6HF + 6H'^{Fo - F) ' 

(14) 



Integrating Eq. ([8]) gives 



Pm^iFonil';>H^{l + zf 



(15) 



where z = ap/a — 1 is a redshift and fji? is a present 
energy fraction of the non-relativistic matter. Then by 
using Eqs. ^ and PU)) the equation of state wue can be 
expressed as 



WDE 



3r- (H-z)(dr/dz) 
3r~3dm\l + zf 



(16) 



where F = df/dR, G^^ is an Einstein tensor, and T^™ is 
an energy- momentum tensor of the pressureless matter. 

In a flat Fricdmann-Robertson- Walker (FRW) metric 
with a scale factor a, we obtain the following background 



where r = H'^{z)/H^. This is the same relation as the 
one derived in Einstein gravity Thus wde is con- 
strained in the usual way from SNIa observations. From 
Eq. p4p we find that the evolution of iode depends upon 
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the models of dark energy. Hence one can test the viabil- 
ity of the models by confronting wde with observations. 

If the scalar field (j) is minimally coupled gravity, e.g., 
/ = R/8ttG + 2p((/), X), the structure of the Lagrangian 
density p{4>, X) can be reconstructed by the evolution of 
the Hubble parameter H{z) [2^. For the models where 
the field (/) is coupled to gravity or the models in which 
the Lagrangian includes non-linear terms in i?, we need 
additional information to determine the strength of grav- 
itational couplings. This can be provided by considering 
the evolution of matter density perturbations. 

III. PERTURBATION EQUATIONS 

We consider the following perturbed metric with scalar 
metric perturbations $ and in a longitudinal gauge: 

ds^ = -{l + 2<^)df +a^{l-2^')d^jdx'dx^ . (17) 

I 



where A: is a comoving wavenumbcr. 

We define the gauge-invariant matter density pertur- 
bation Sjn, as 

S„i = ^™ + 3Hv , where v = avm . (23) 

Pm 

Then Eqs. (Ell) and 1^ yield 

Sm = ^v + 3{^ + Hvy, (24) 

u = $ , (25) 
from which wc obtain 

Srn + 2HS„, + ^$ = W + 6HB , (26) 

where B = + Hv. 

Following the approach in Rcf . [l^, [13 , we use a sub- 
horizon approximation under which the leading terms 



We decompose the field into the background and inho- 
mogeneous parts: (/> = (f){t) + 5(/)(i,x). In what follows 
we drop the tilde for simplicity. The energy momentum 
tensors of the non-relativistic matter are decomposed as 
^0* = -(Pm + Sp„i) and T° = -pmVm.a, whcrc Vrn IS a 
velocity potential. 



The Fourier transformed perturbation equations are 
given by [2^ 



(18) 

(19) 

(20) 
(21) 
(22) 



correspond to those containing and 5m (or Spm) in 
Eq. ([^ and also in Eqs. ((l8l)-([T9l). Basically the terms 
on the r.h.s. of Eq. ([^S)) give the contribution of the order 
H^'^ , which implies that they are negligible relative to 
the term (fc^/a^)<I> for the modes deep inside the Hubble 
radius (fc^ > a?H'^). 

If the mass of the field perturbation (50 is larger 
than the term k/ a then we need to take into account this 
mass term. The expression of was derived in Ref. [20| 
in coupled dark energy models with the Lagrangian den- 
sity p{(j}^X). In Einstein gravity with a standard scalar 
field the mass squared is given by to^ = V^^^ — 2(j)'^. 
When the field (j> is responsible for dark energy the terms 
Vcficl, and 2(p^ are of order or less, which then gives 
i$ H. Hence the approximation neglecting the mass 
term relative to k/a is justified in such a model. 
There may be some specific f{R,(j),X) models in which 
the condition \m^\ ^ k/a is violated, but we do not con- 



3iJ(* + i/$) + ^* + -3- 
21" 



-\{f,^H + f.xSX) + ^4,^{lx<pS(b + f,xxSX) + lx<pS<P - 3HdF 



3H^ + 3H--^]6F + 3F{^ + H^) + {3HF - f,xr)'i> + Sp„ 



= 0, 



f V \ • k'^ 

3H + ^]S(I)+ —5<i) - 0(3^' + $) 



X 



- 2/,^$ + ^{a^5f^xy ~ Sf,^ = , 



* = $ + — 



5pm + iHSpm = Pm 3^* 



a 
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sider such cases. 

Then Eq. (j26p is approximately given by 



fc2 



0. 



The next step is to express $ in terms of i5; 
Eq. (dH) wc find 



Ehminating the term 5* by using Eq. (PD)) gives 



(27) 
„. From 



fc2 1 
'2^~F ~ 2F 



Sp„ 



(28) 



(29) 



In what follows we shall study the case in which F 
depends on </> and R but not on X , i.e., 



F^F{cp,R). 



(30) 



This actually includes most of dark energy models pro- 
posed in the current literature. Then 5F in Eq. p9|) is 
given by 



SF = F 



FrSR, 



(31) 



where SR is 
SR = 2 



fc2 .\ pi 
-3H] $-2—5' 



(32) 



-2— f $ + 2— 



Again we used the fact that the first five terms in Eq. 
are of order H^'^ or less. Plugging Eq. (|3ip into 

Eq. dsn), we find 



2fc2 $. 



2F„ 



(33) 



Taking notice that 5R includes the term k'^/a^, the 
quantity 6f is approximately given by 



Sf = l^H + f,x5X + f^R8R 
~ F5R. 

Then from Eq. p9p we find 
a'' 

which leads to the following relation 



(34) 



(35) 



(36) 



Plugging this into Eq. ([32|) . we get 



6R' 



.X 



fe2 F.H 



(37) 



From Eqs. §^ and JST]) the gravitational po- 

tential $ is expressed as 



/,x+4(/,x|^^ + T^ 



2F 



/.x+3(/.x^^ + ^ 



(38) 



where we used — 5pm/pm under the sub- horizon ap- 
proximation. Hence the equation (|27p of matter pertur- 
bations yields 



2H6,n ~ AnGcSPnAn ^ , 



(39) 



where the effective gravitational "constant" on sub- 
horizon scales is given by 



G, 



off 



1 lx+Mlx^^ + ^ 



BttF 



/.x + 3(/,xi^ + ^) 



(40) 



From Eq. p8|) the gravitational potential is 



(41) 



which corresponds to a Poisson equation in the Fourier 
space. In what follows we use the standard equality 
rather than the approximate equality (~) for the results 
obtained under the sub-horizon approximation. 

We also define a parameter rj that characterizes the 
strength of an anisotropic stress: 



V 



$ - * 



Using Eqs. ^ and ^ we obtain 



2/, 



X 



F 



' " J ^ ~~F J ' F 

The gravitational potential ^' satisfies 
/x + 2 (/.xfr-^ - 



(42) 



(43) 



2F 



f.x + 3{lx^^ + ^ 



(44) 



We define another parameter q via (k^/a^)^ ~ 
—iTrGoqpmSrm where Go is a gravitational constant mea- 
sured in the solar system experiments today. Then q is 
given by 



1 /.x + 2 /.xf 



" J,X + 'i[^J,X-^— + — 



(45) 
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Defining a combination of parameters, S 
we obtain 

1 



9(1 + '7/2), 



(46) 



This agrees with the result in Ref. [l^ derived in the 
specific scalar-tensor model: / — F{(p)R + 2X — 2V{(f))^ . 
For this model the parameter 77 reduces to 77 = F'^/ {F + 
i^^), which again agrees with the result given in Ref. [l3l |. 

In order to confront the modified gravity models with 
the observations of weak lensing, we use the fact that the 
potential that characterizes the deviation of light rays 
corresponds to $wL = $ + 4" From Eqs. ((SHI), (011) 

and (|46p we find that the lensing potential satisfies 



4> 



WL 



iT^Go—pmSm^ ■ 



(47) 



The effect of modified gravity theories manifests them- 
selves in weak lensing in at least two ways. One is the 
multiplication of the term E on the r.h.s. of Eq. P7)) . 
Another is the modification of the evolution of 6m due 
to the change of the effective gravitational constant Goff . 
The growth index 7 of matter perturbations is linked to 
the parameters E and 77 p^ . Thus two parameters (E, ?/) 
will be useful to detect the signature of modified gravity 
theories from the future survey of weak lensing. 

It can happen that the scales of weak lensing are in the 
region of non-linear clustering, in which case we need to 
map the linear power spectrum of the lensing potential 
into a non-linear one. In the context of modified gravity 
theories mapping formulas have not been well known. 
We leave the analysis of such non-linear regimes in weak 
lensing for future work. 

In Einstein gravity with the Lagrangian density / = 
R/8ttG + 2p{<j>,X) we obtain the standard equation for 
matter perturbations: 



0, 



(48) 



where a prime represents a derivative with respect to 
= In a, and 

If WcS and flm are constants then the solution for Eq. ((48l) 
is 



One has w^s — and rj,„ — 1 during a matter-dominated 
epoch provided that the contribution of the scalar field is 
negligible. Hence the matter perturbation grows as 6m cx 
a cx t^^^. However the evolution of 6m is modified once 
the energy density of the scalar field becomes important 
relative to the matter density. 

At the end of this section we consider the Brans-Dicke 
theory [2|: 



(52) 



where lobd is a Brans-Dicke parameter. In this case the 
effective gravitational constant is given by 



1 4 + 2lubb 
(j)3 + 2wbd 



(53) 



which agrees with the result in Ref. |30[ . We have Gcff 
1/^ in the General Relativity (GR) limit (wbd 00). 
The deviation from GR is significant when wbd is not 
much larger than unity. When ojbd = 0, for example, 
the effective gravitational constant is 4/3 times larger 
than that in the GR case. This modifies the evolution of 
matter perturbations. However local gravity experiments 
place the bound on the present value of the Brans-Dicke 
parameter as cjbd.o > 4x 10'' [Slj. This shows that unless 
ciJbd is very much smaller than the present value during 
the matter epoch it is difficult to see the signature of 
modified gravity in the large-scale structure formation. 



IV. f{R) GRAVITY 

In this section we study the evolution of matter per- 
turbations in modified gravity theories where / is the 
function of R only. In this case the effective gravitational 
constant is given by 



1 1 



SnFl- 



where 



RFj 



(54) 



(55) 



6m = C+CL^^ + c_a" 



(50) 



where 



"± = 4 



iw,s - 1 + ^{iw,s ~ 1)2 + 24^!™ . (51) 



^ In Ref. rOf the authors used a dimensionless function F = SttGF, 
where G is a bare gravitational constant. Then one has S = 
G/FGo. 



The parameter m was first introduced in Ref. [2J]. 
This characterizes the deviation from the ACDM model 
(/(i?) = R/8ttG — A). The anisotropic parameter 77 is 
given by 



(56) 



In what follows we shall consider two different situa- 
fc^ ™ ^ 1 fc^ 



tions: (i) -^m ^ 1 and (ii) 



771 < 1. 
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Tin > 1 



^ , which thus corre- 



In this case one has G^ff '— o— e?^ 
sponds to Brans-Dicke theory with o^bd = [33[, see 
Eq. ([55]) . Note that the anisotropic parameter is of or- 
der unity (77 ~ 1) in such a case. Since the condition 

/a^R 3> 1 holds under the sub- horizon approximation 
{k 3> aH), one can in fact reahze -^r^m ^ 1 provided 
that TO is not very much smaUer than unity. Then from 
Eq. (|48p the matter perturbation equation is approxi- 
mately given by 



0^ 



5'r„ - 2f7m5„ 



0, 



(57) 



where and are defined in Eq. 

In order to estimate the evolution of 5m and $ analyt- 
ically, let us consider the constant to model, i.e., 



f{R)^aR'+"'-A, 



(58) 



where a and A are constants. In Ref. |24l | it was shown 
that the matter-dominated epoch corresponds to a fixed 
point "P5" satisfying 

m o -, to(7+ IOto) 

Wcff^-— , rJ™ = l-— — r^, 59) 

1 + TO 2(1 -|- to)^ 

where |to| <C 1. 

Plugging Eq. ([55]) into Eq. ([57)) . we obtain the solution 
for Sm in the form (l50l) with 



-(1 + 4m) ± V(3 + 4to)(11 - 28^ ...A 

'^i = T7T— — X • (60) 

4(1 -I- to) 

Since the growing mode corresponds to the power-law 
index n+, the matter perturbation evolves as 



where 



^(3 + 4to)(11-28to) - 1 - 4to 
6 



(61) 



In Ref. [2J| it was found that the viable matter epoch 
exists only for positive to close to 0. In the case of neg- 
ative TO the matter point P5 is unstable against pertur- 
bations around the fixed point. When to = —1/4 one 
has uicff = 1/3 and Q„i = 2. which corresponds to a 
matter-dominated epoch (0MDE) [s^. For the models 
f{R) = R- l3/R^ {n > 0) it was shown in Ref. [H that 
the standard matter era is replaced by the ^MDE. From 
Eq. (|60p we find that the matter perturbation evolves as 
6m oa oa t, which grows more rapidly than in the stan- 
dard case {6m oc a cx t^/'^). 

The power-law indices n+ and ri+ are positive for 
< TO < (\/73 — 3)/16, whereas they are negative for 
(773 - 3)/16 < TO < 11/28. When to > 11/28 the mat- 
ter perturbation exhibits a damped oscillation. Both n+ 



and n+ get larger as m decreases to zero. However we 
have to caution that we can not take the limit to ^ 



because of the breakdown of the condition 



TO > 1. 



In this limit the evolution of 6m is no longer described 
by the solution ((6T|) . 

From Eq. (|41]) we find that the gravitational potential 
evolves as 



$oct-, p,^^(^+^-)(^^-/^-)+^-'^(63) 

One has p+ > for < to < 1/4 and p+ < 
for 1/4 < TO < 11/28. Hence $ is not constant ex- 
cept for the special case to = 1/4 and the 0MDE case 
TO = —1/4. The variation of the gravitational potential 
leads to an Integrated-Sachs- Wolfe effect in the CMB 
spectrum. Thus it should be possible to constrain the 
magnitude of to from CMB observations. 



B. 



Let us next consider the case in which the condi- 
tion, -A-^m ^ 1, is satisfied on the scales around which 
large-scale structure is formed. In this case one has 

Gcff — ^1 -t- -^^'m^ and 77 ~ The matter 

perturbation equation is approximately given by 

fc2 



1'^ 



a?R 



0. 



(64) 



Since the condition to ^ -titto <C 1 holds under the sub- 
horizon approximation, we only pick up the correction 
terms that contain -^to. 

Expressing the solutions of this equation in the form 
6m = exp(J ujdN) and using the approximation <C 



we obtain the growing-mode solution 

3 k^ 
"+ = ^+5^' 



(65) 



If TO is constant, the second term on the r.h.s. of this 



(62) equation is proportional to a 



during the matter 



era. Then the evolution of the matter perturbation is 
given by 



oc a 



1+1/3 



OC t3 



(l+f/3) 



where 



^— — TO = C/C TO- , 

a^RJSl ma 



(66) 



(67) 



Here we have introduced a constant C satisfying the rela- 
tion l/(a^i?) = Ce^ . The gravitational potential evolves 
as 



4>oc {\ + (3)tif^ . 



(68) 



In the limit /3 ^ one obtains the standard result: 
6m t^/^ and <I> = constant. When /3 is positive, the 
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growth rates of Sm and $m are larger than in the stan- 
dard case. If (3 grows to the order of unity, the results 
([66|l and ((68|) arc no longer valid. 
In order to satisfy the local gravity constraint (LGC), 

7 2 

we require that the condition -^r^m ^ 1 holds at the 
present epoch. The severest constraint may be obtained 
by laboratory experiments in which a strong modifi- 
cation of gravity is not observed on the scales up to 
Afe ^ a/fc ~ 1 mm. This gives the following constraint 



m(z = 0)<(Afc/Fo"i)2 



10 



-58 



(69) 



SCALAR-TENSOR GRAVITY 



One can study the evolution of perturbations in the 
Jordan frame as we have done in /(i?) gravity models. 
Alternatively the dynamics of perturbations can be un- 
derstood by making a conformal transformation to the 
Einstein frame. This is particularly useful when we re- 
late scalar-tensor models with coupled DE scenarios [llj 
extensively studied by many authors (see [6| for refer- 
ences). We shall make a conformal transformation jSTj 



(73) 



S 



where we used R-^ and H^^ - lO^^ cm. Note that 
this agrees with the result in Ref . [l^ that was derived by 
using the effective mass of a scalar-field potential in the 
Einstein frame (sec Rcfs. [lH for recent works of LGC 
in f{R) gravity). The condition ([55)1 can be weakened 
by taking into account the fact that the scalar curvature 
R in the regime of a local structure such as an earth is 
much larger than the cosmological one. In solar system 
experiments the scale corresponds to a value around 
Xk ~ 1 au. which is much larger than in the case of lab- 
oratory experiments. In such cases the constraint on m 
becomes much weaker than the one given in Eq. (|69p . 
although it is not easy to obtain the values of m close to 
the order of unity [3J| . 

The variable m generally changes with time apart from 
the model / — aR^'^"^ — A discussed above. One can 

.2 

consider models in which m satisfies ■^rji'm 3> 1 during 

,2 

the matter epoch and then enters the regime -^r^'m ^ 1 
in the dark energy era with the decrease of m. The 

.2 

anisotropic stress parameter rj decreases from 1 to 2^^m where 
together with the change of the quantity E given in 
Eq. It will be of interest to place observational 

constraints on such models by using the future data of 
weak lensing as well as the LSS data. 



where a tilde denotes quantities in the Einstein frame. 
Then the action in the Einstein frame is given by 



1 
2k 



2 i? + ]5(0, A) +£„(</.) 



(74) 



where = SttG (G is a bare gravitational constant) and 
p(^, X) = U^\' X + Pi^, FX) . (75) 



In what follows we shall use the unit 



= 1, but we 



restore the gravitational constant G when it is needed. 
We also have the following relations 



= P^/F\ (76) 
H) is transformed as 



5 = VFa, dt = VFdt, 
Then the continuity equation ([ 



dt 



dt 



(77) 



(78) 



Hence the non-relativistic matter is coupled to the field 
(j) through the coupling Q{(j)). 

The perturbations in the Einstein frame are related to 
those in the Jordan frame via [2a| 



In this section we shall consider scalar-tensor gravity 
with a Lagrangian density 



/(i?,0,A) = F(0)i? + 2p(0,A). 



(70) 



Note that this includes most of scalar-field dark en- 
ergy models such as quintessence 0, k-essence [1| and 
tachyons [s^ . In this case the effective gravitational con- 
stant Geff and the anisotropic parameter 77 are given by 



V 



1 2p,x+AFl/F 



SttF 2p,x + iF^/F 
Fl 



P,xF + Fi 



(71) 
(72) 



In the case of quintessence with the Lagrangian density 
p = X — V{(j)) the above results agree with those obtained 
by Boisseau et al. [131 and by Amendola et al. [isj . 



25F 



$ = $ 



26F 



(79) 



Under the sub-horizon approximation one can regard 
(5m — 5rm whereas the 2SF/F contribution can not be 
neglected for the gravitational potential. Using the rela- 
tion P X ^ ("^) ~'~ ^ together with the definition of 
Q given in Eq. (|78p . we find that Eq. ([M)) is written in 
terms of the quantities in the Einstein frame: 



^L + i/(2- 



^ H dtjdt 



-AttG 1 



-r PmOm = . 



m 



In this section we use a prime to represent a derivative 
with respect to the number of e-foldings N = ^ H dt 



8 



in the Einstein frame. Then the above equation can be 
written as 



to the case g{Y) = 1 — ci/Y. For the models 
have the foUowing relations during the (pMDF, [40 



we 



2 + -^ 
H 



2Q^ 

P.x 



(81) 



where Clm = SnGprn/SH^ ■ This fully agrees with the re- 
sult in coupled dark energy scenarios derived by Amen- 
dola [l^l without any reference to the Jordan frame (see 
also Ref. (H). 

Equation (|80p shows that the effective gravitational 
constant in the Einstein frame is given by 



Gcff = G 1 




2px+ ^F\/F 



From Eqs. and l|79|) the gravitational potential in 
the Einstein frame is 



2p.x + 3FyF 



2p^x + 4Fl/F ' 
which satisfies the relation 



(84) 



The effective gravitational potential acting on the matter 
is not $ but l>* = $ + Q5cj) = $ H^], i.e., that in the 
Jordan frame. In fact, from Eqs. ([5^ and ((51)). we obtain 



-ATrGcff-r^PmSn 



(85) 



In order to see the effect of an interaction between the 
field (j) and the matter analytically it is convenient to 
study the constant Q case, i.e.. 



F(0) 



-2Q4> 



(86) 



In the case of an ordinary field with an exponential poten- 
tial [p{(f>, X) = X — ce~'*"^], it is known that there exists a 
0MDE scaling solution satisfying fl^ = Wcs = 2Q^/3 = 
constant [2lj . More generally the existence of scaling so- 
lutions restricts the form of the Lagrangian density to 
be 



(87) 



where a is an arbitrary function and A is a constant quan- 
tity [la. It was further shown in Ref. ^ that the 0MDE 
exists for the models of the type 



giY) ^co + J2 <^nY- 



(88) 



?i>0 



where cq (> 0) and c„ are constants. Note that 
quintessence with an exponential potential corresponds 



2^ 

3co ' 



2Q . 
' P X 

Co 



Co . 



(89) 



where Sl^ is an energy fraction of the scalar field satisfy- 
ing -I- ri„i = 1. 

Then from Eq. (|8ip . the solution for matter perturba- 
tions during the 0MDE is given by the form ([50]) with 
the power-law indices 



= 1 



2g^ 

Co 



71_ = 

2 



Co 



(90) 



Since the scale factor grows as a oc 1 3=0+20^ ^ the evolution 
of the matter perturbation in the Einstein frame is given 
by 



1 20^ _2co+4Q^ 

5m OC a "^0 (X t3co+2Q'2 



(91) 



(83) From Eq. jHS]) we obtain 



4Q 



Int. 



^ 3co + 202 

Using the relations (|76p between two frames we find 



(92) 



t (X i3co + 2Q2 



(93) 



Hence the evolution of matter perturbations in the Jor- 
dan frame is given by 



Sr,i OC a^-o-iQ'-' OC <3co-2Q^ 



(94) 



Thus, in the presence of the coupling Q, the growth rate 
of matter perturbations is larger than in the case of Ein- 
stein gravity. From Eqs. ([M]) and ([SS]) we find that the 
gravitational potential is constant in both Jordan and 
Einstein frames: 



$ cx i' 







^cxP, 



(95) 



which is a rather peculiar property of the (pMDF,. Recall 
that this property also holds for the 0MDE solution in 
the /(i?) gravity. 

From Eq. (j75|) the Lagrangian density in the Jordan 
frame corresponding to the 0MDE solution is given by 



R + 2{co-6Q^)X 



^2cnX^- 



ng-n(A+2Q)</. 



n>0 



(96) 



Thus the field (p has a universal coupling e"^'^'^. For the 
model g{Y) ~ l — ci/Y the above Lagrangian density can 
be viewed as the dilaton gravity with an exponential po- 
tential y((/)) = 2cie"(^+2'5)<^. It is interesting that string 



9 



theory can give rise to the 0MDE solution along which 
f2m and Weff are constants in both Jordan and Einstein 
frames. We note that Q is required to be smaller than the 
order of unity to reproduce a standard matter era [2]| . 
whereas string theory typically provides the coupling Q 
of order one at the perturbative regime [4l|. In the re- 
gion 0^1, however, the coupling may become weak as 
in the context of a runaway dilaton scenario (4^ . 

VI. CONCLUSIONS 

In this paper we derived the matter perturbation equa- 
tion ((39)) with the effective gravitational constant ([40]) for 
a Lagrangian density f{R, <j), X) without a direct cou- 
pling between R and X. This analysis covers most of 
modified gravity models proposed in the current liter- 
ature and will be useful to detect the deviation from 
the ACDM model from the future surveys such as weak 
lensing and LSS. We have also evaluated the anisotropic 
parameter rj and the quantity E = q(l + ri/2) in order 
to confront the models with future observations of weak 
lensing. see Eqs. (P5|) and (|1S)) . We have applied our 
results to (i) f{R) gravity and (ii) scalar-tensor gravity 
with the Lagrangian density f{R) = F{4))R + 2p{(j), X). 

In f{R) gravity models the effective gravitational con- 

.2 

stant has a scale-dependent term -^^Tn, where m — 
Rf ,rrI f .R characterizes the deviation of the ACDM 
model. The local gravity constraint is satisfied for the 

,2 

models in which the condition, -r^m <C 1, holds for 
the scale of the order a/fc = 1 mm at the present epoch. 
If we take the cosmological value R ~ H^, this gives 
a very stringent constraint: m{z = 0) ^ 10~^®. This 
is weakened by using a local value of R much larger 
than Hq. One can consider models in which the con- 

,2 

dition, -^r^m ^ 1, holds during the matter epoch on 
the scales around which large-scale structure is formed. 
For the constant m model (/ = Q;i?^+™ — A) we have 
analytically derived the evolution of matter perturba- 
tions 5m as well as gravitational potentials <& during the 
matter-dominated epoch. Even when to <IC 1 this is dif- 
ferent from the evolution in Einstein gravity {5m oc a and 
$ ^constant), which will be useful to place constraints 
on the value m from future high-precision observations. 

The scalar-tensor gravity with the Lagrangian density 
f[R) = F{(j))R + 2p{(j), X) correspond to coupled DE 
models in the Einstein frame with a coupling Q{(j)) = 
—F^^/2F between the scalar field and dark matter. We 
reproduced the equation of matter perturbations in cou- 
pled DE scenarios with the k-essencc Lagrangian den- 
sity p{(p, X) by making a conformal transformation to 
the Einstein frame. Since the evolution of perturbations 
in coupled DE models has been extensively studied in lit- 
erature, it is convenient to pay attention to the relation 
between Jordan and Einstein frames in order to discuss 
perturbations in the scalar-tensor gravity. In fact, for the 
models in which the matter dominated epoch (0MDE) 



is present, we analytically derived growth rates of per- 
turbations in both Jordan and Einstein frames. We also 
obtained the form of the Lagrangian density in the Jor- 
dan frame giving rise to the (/)MDE solution. 

The difference between f{R) gravity and scalar-tensor 
gravity may be understood in the following way. Taking 
into account a mass M of the perturbation in the field 
(j), the effective gravitational constant in Eq. ([5^ is given 
by M 

Geff=.G(^l + ge~^«^ , (97) 

where £ is a length scale. Here wc assumed that the sound 
speed Cs of the field 4> is of order unity. In scalar-tensor 
models, even if the mass M is very light as M ^ Hq, the 
second term on the r.h.s. of (P7|) can be much smaller 
than unity to satisfy the local gravity constraints by 
choosing a small coupling = {F^/2F)^ ^ 1 (pro- 
vided that p is of order one). In f{R) modified gravity 

models, however, the coupHng Q is fixed as Q = —I/VQ 
[33 |. Hence we have to choose a heavy mass ~ ^/P,R 
to satisfy the local gravity constraints, i.e., e~^^^ ^ 1 
for the scale £ ~ 1 mm [34| . In fact this is equivalent 
to choosing very small values of to given in Eq. (j69p . 
In scalar-tensor gravity the coupling-dependent term is 
more important than the scale-dependent term provided 
that AI <C fc/a, whereas in f{R) gravity the scale- 
dependent term plays a crucial role because of the fixed 
large coupling Q. 

Although the perturbation equation we derived can be 
applied to many modified gravity models, it does not 
cover the DE models in which higher-order curvature cor- 
rections such as a Gauss-Bonnet (GB) term are present 
[ist . In Ref. [11] the equation of matter perturbations 
was derived in the presence of the GB term coupled to 
a scalar field (f> to p lace constraints on GB DE models. 
A recent paper [43 shows that the GB energy fraction 
in the present universe is severely constrained by solar 
system tests. Moreover it is known that tensor perturba- 
tions typically show negative instabilities if the GB term 
is responsible for the accelerated expansion [i^ . It will be 
of interest to extend our analysis to more general models 
that include such higher-order curvature corrections. 
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